Physical Limits of the ballistic and non-ballistic Spin-Field-Effect Transistor: Spin 

Dynamics in Remote Doped Structures 
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We investigate the spin dynamics and relaxation in remotely-doped two dimensional electron 
systems where the dopants lead to random fluctuations of the Rashba spin-orbit coupling. Due 
to the resulting random spin precession, the spin relaxation time is limited by the strength and 
spatial scale of the random contribution to the spin-orbit coupling. We concentrate on the role 
of the randomness for two systems where the direction of the spin-orbit field does not depend 
on the electron momentum: the spin field-effect transistor with balanced Rashba and Dresselhaus 
couplings and the (Oil) quantum well. Both of these systems are considered as promising for the 
spintronics applications because of the suppression of the Dyakonov-Perel' mechanism there makes 
the realization of a spin field effect transistor in the diffusive regime possible. We demonstrate 
that the spin relaxation through the randomness of spin-orbit coupling imposes important physical 
limitations on the operational properties of these devices. 



I. INTRODUCTION 

The simultaneous control of the spin and charge dy- 
namics in low-dimensional semiconductors by electrical 
means has been one of the major goals of the semicon- 
ductor spintronics field for over a decadei^ A strong 
focus within this field has been systems with spin-orbit 
(SO) coupling, which connects both the spin and charge 
degrees of freedom. In strongly asymmetric heterostruc- 
tures grown along the [001] direction, the dominant SO 
coupling is of the Rashba type due to the inversion asym- 
metry of the confining potential;^ and has the form 
Hr = cxn(a x py — o~ y p x )/h, acting as a momentum depen- 
dent Zeeman field that aligns the spins of electrons per- 
pendicular to their in-plane momentum. Here, a are the 
Pauli matrices, and p = (p x ,p y ) is the in-plane momen- 
tum. The controlled precession of injected spin polarized 
carriers around this effective Zeeman field produced by 
the Rashba coupling in momentum space is the basis for 
the seminal proposal of Datta and Das of a ballistic spin- 
field-effect transistor (DD-SFET) 5 . However, a decade 
long quest to realize such a device has yielded disappoint- 
ing resultsi^ In part, this is due to the fact that the 
SO coupling, which crucially influences spin transport, 
also determines the spin relaxation rate in weakly dis- 
ordered systems through random precession due to elas- 
tic scattering processes, the so called Dyakonov-Perel' 
(DP) mechanism^ hence limiting the validity of the two- 
dimensional DD-SFET to the ballistic regime. 

Schliemann et ali have recently suggested a way to de- 
sign a SFET in (OOl)-plane based structures that should 
work on the diffusive regime, hence without the need to 
reach this ellusive ballistic regime, by using the interplay 
between the Rashba and the Dresselhaus SO coupling, 
Ho = aoicyPy ~ cr x p x )/H, which originates from the 
lack of inversion symmetry in the unit cell of the host 
material^ The key aspect of this proposed device is 
that the momentum dependence of the direction of the 



effective Zeeman field is removed if the Rashba and Dres- 
selhaus coupling strengths are equal with cxr = ±otr>. In 
this case the in-plane SO field is parallel either to the 
[110] or [110] direction and does not depend on the direc- 
tion of the electron momentum. It was suggested recently 
m Refs-HHU that the freezing of one of the electron mo- 
mentum components in one-dimensional channels sup- 
presses the DP mechanism of spin relaxation and leads 
to a much longer spin lifetime than in the 2D quantum 
wells. This idea, is, however, different from the proposal 
of Ref . , where the DP mechanism becomes much less 
efficient through the combined effects of the Rashba and 
Dresselhaus terms (see also in Ref. [l2| ) . 

By decoupling the momentum and the spin- 
quantization axis orientation, the DP mechanism is 
greatly reduced, and hence the spin coherence time in- 
creased making such a device apparently realizable in 
both the diffusive and ballistic regimes. The condition 
ac = an is indeed a possible one since it has been 
recently demonstrated that by applying external bias 
across the quantum well, one can change the magnitude 
of a R in InGaAs/InAlAs-basediii^ and GaAs/AlAs- 
basedl&J&JZ„ structures and reach ratios of ao/ctR close 
to one. Another example of the decoupled SO field and 
momentum directions is provided by (Oil) zincblende 
quantum wells (QW), where the Dresselhaus SO coupling 
leads to the effective field perpendicular to the well plane 
at arbitrary electron momentum, as discussed later. 

Here we consider an aspect of the practical realization 
of this device which will have important limiting conse- 
quences in its operating properties. We show that the 
randomness of SO coupling arising in doped systems de- 
stroys the regular precession of electron spins and leads to 
a spin relaxation both in the SFET device and the (Oil) 
quantum wells even in the ballistic regime. We show that 
the size of the devices based on the idea of this effective 
SO decoupling is limited even in the best possible real- 
izations of the devices and, therefore, the spins cannot be 
coherently transferred at arbitrary distances. In most of 
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the low-dimensional QW-based semiconducting systems, 
the SO coupling and finite concentration of electrons are 
achieved either by one-side or asymmetric two-side re- 
mote doping at distances between the symmetry plane of 
the QW and the dopant layer of the order of 10 nmiiiiii^ 
The spatial non-uniform distribution of dopant ions leads 
to fluctuations of the SO coupling, and therefore, the con- 
dition of cud = cur in a SFET can only be achieved on the 
average throughout the sample; the fluctuations of this 
coupling are analogous to the classical "shot" noise in 
the transport properties. The effect is, in some aspects, 
similar to that investigated by Mel'nikov and Rashba in 
doped bulk Si and Ge— At the same time, recent the- 
oretical analysis^ shows that in the case of strong SO 
coupling (large an) achieved in highly asymmetric het- 
erosructures the condition old — ur is hardly achiev- 
able by applying a bias, and therefore, we are left with 
the quantum wells where the moderate Rashba coupling 
arises due to the doping or weak structural asymmetry. 
For a comparison to possible experimental realizations, 
we shall consider InGaAs structures, where the Rashba 
coupling can be reliably influenced in a wide range by 
the applied bias, as has been proven in Refsiiiii 

We organize the rest of the article as follows: In Sec. II 
we discuss the properties of the SO coupling considered. 
In Sec. Ill, we discuss the spin relaxation in the proposed 
field-effect (FET) with balanced Rashba and Dressclhaus 
SO coupling terms and in Sec. IV, the one for the (Oil) 
quantum well. In Sec. V, we present conclusions for these 
results. 



where e is the dielectric constant, and the summation is 
performed over the dopant layer. 

Fig 1. presents the pattern of an(p) at the spatial scale 
of 200 nm for zq — 20 nm obtained by a Monte-Carlo 
produced "white- noise" distribution of the dopant ions^. 
The fluctuations of an(p) correlated on the spatial scale 
of the order of zq become smaller and smoother with the 
increase of the distance to the dopant layer, as is shown 
by one-dimensional cuts in Fig. 2. 

The Rashba parameter is a sum of the mean value 
produced by the asymmetric field (a) and a random con- 
tribution with zero mean: 

a R (p)={a)+a I (p), (2) 

with (a) = ao + 27raso (e 2 / 6 ) ^> where n is the 
mean value of the 2D dopant concentration. For the 
S(z + zq)— doping in quantum wells with a$ = the 
relative amplitude of the fluctuations (a 2 ) / (a) — 
1 / y/ 87mzg ,22*24 In contrast to the regular term, which 
demonstrates a dependence or(V) on the applied bias 
V across the well, this random contribution remains rel- 
atively weakly influenced and, therefore, cannot be re- 
moved. 



II. RANDOMNESS OF THE ZEEMAN SO 
FIELD. 

We consider a two-dimensional (2D) channel and a 
very thin dopant layer (<5-doping) separated by a dis- 
tance zq with a two-dimensional concentration n(r) = 
^2jS(r — Yj) of dopants at points r, with charge e and 
r = (r x , r y ) being the 2D in- plane vector. 

Since we are interested in quantum wells, where the 
Rashba SO coupling arises due to an asymmetric dop- 
ing, we assume that the effect of the electric field of 
dopants on the SO strength is linear in the z— component 
of the field E z (p), that is an{p) = asoeE z (p). Here p is 
a point at the plane of a sufficiently narrow 2D chan- 
nel and «so is a phenomenological system-dependent 
parameter—. This relationship takes into account the 
linear in E z polarization of the electron wavefunction 
as well, as found in Ref. [2l|. In more asymmet- 
ric heterostructures or asymmetrically grown quantum 
wells, the structural asymmetry will lead to a regular 
coordinate-independent shift ao to the Rashba term, 
which has no influence on the effects of the randomness. 
The z-component of the Coulomb field due to the dopant 
ions is given by 

E z (p) — —— 77 oTTTTT, (1) 

e Y [(p-r,) 2 + z 2 ]3/2' 



III. SPIN RELAXATION IN THE 
NON-BALLISTIC SFET CONDITION. 

The randomness in the direction of the spin precession 
axis can be removed when at some applied bias V, the 
condition old = ±an(V) is satisfied and a R (V) is spa- 
tially uniform. Therefore, the direction of the linear in 
the in-plane momentum Zeeman field does not depend 
on the momentum directions, and the spin precession 
becomes more regular. Here a new integral of motion 
£ = (cr x + <7y) j\/2 appears and supports the regular- 
ity of spin dynamics i We emphasize here that the exact 
matching is possible only at the limit of zero electron mo- 
mentum, where both the Rashba and Dresselhaus Hamil- 
tonians are linear in p. At finite momenta the Dressel- 
haus Hamiltonian must be supplemented by the bulk fc 3 - 
like term, which in the (001) quantum well has the form: 
a c k y k x (a x ky — cr y k x ), where a c is the Dresselhaus con- 
stant for the bulk, and k = p/h is the wavevector. The 
role of this term depends considerably on the quantum 
well width. This effect, as well as the possible but not 
yet well known momentum dependence of the Rashba pa- 
rameter, can hamper the ability of design the SFET even 
if the linear contributions are balanced. 

The total Hamiltonian describing the system written 
in the Hcrmitian form: 
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Hso 


— Hd + Hr, 


+ \ Wx {k y , {a R (p) - {a R })} - G y {k x , (a R (p) - (or))} 


H R 


= (aR)(cr x ky - o y k x ) 


H D 


= ao{(Jyky - a x k x ) -i 





where {fcj,aij(p)} stands for an anticommutator. The 
Dresselhaus coupling constant ol-q — a c (k^), where (k 2 ) 
is the expectation value of the k 2 = —d 2 /dz 2 operator in 
the ground state, and a c « 25 eVA 3 in GaAs and InAs- 
based structures 2 ^*. In the case of the rigid walls of the 
QW with the width w, {k 2 z ) = (tt/w) 2 . 

Both the cubic terms in the Hamiltonian and the ran- 
domness of SO coupling lead to the spin relaxation in 
the SFET device. Our main interest in this paper is 
the role of the randomness. However, to find the regime 
where the effect of the randomness dominates, we will 
first discuss the contribution of the a c k y k x (a x k y — a y k x ) 
term for the spin relaxation in the quantum well in the 
case an — (ap) . We will use the approach developed in 
Refs. 6,8 '— , where spin relaxation is described by: 

3 

with Tij being the anisotropic spin relaxation tensor. For 
the relaxation of E (expected to be conserved in the an — 
(ap) case) we obtain: 



where the contribution of the k 3 — originated terms yields 
Td — a 2 k e Tk/8H 2 , with being the momentum relax- 
ation time. We will demonstrate that under realistic con- 
ditions, this term is much smaller than the contribution 
of randomness of the SO coupling, which we consider 
next. 

The randomness of the SO coupling discussed here 
makes the precession irregular and, therefore, pushes the 
range of parameters in which such a SFET can be real- 
ized toward the ballistic regime. We consider a SFET 
configuration where electrons are injected through per- 
fect nanoscale quantum point contacts^ located at points 
x = 0, y = yi, where i = 1, 2, . . . N numerates the con- 
tacts and the electrons and then move along the x axis 
as shown in Fig.l interacting with the random SO field. 
The initial spin state of all electrons is 

&(i = 0,a: = 0) = -^f * W4 ), (6) 

with all the spins initially polarized along the x + y di- 
rection. The shape of the injected electron density along 
the y-axis, \i^(y — yi)\ 2 formed by the quantum contact 
is conserved during the electron propagation such that 



the corresponding wavefunction has the form: ^i(x, y) — 
tp(y — exp(ik x x). We assume that the y-axis spatial 
distribution of the wavefunction is much less than zo, 
and, therefore, the spin of the zth electron interacts with 
the local Rashba field aR(x,yi) and its spatial deriva- 
tives. To estimate the maximum spin-coherence lifetime 
we assume that the movement is ballistic, that is the 
momentum of injected electron is conserved, and, there- 
fore, only the randomness of aR(x,y) leads to the spin 
decoherence, thus, presenting the spin relaxation mecha- 
nism in relatively high-mobility structures with the elec- 
tron free path much larger than 10 nm. This restriction 
implies that the electron kinetic energy is much larger 
than the potential fluctuations and ignores quantum me- 
chanical interference effects which are small in the larger 
system sizes considered here. Relaxation of this approx- 
imation will further decrease the spin-coherence lifetime 
in the diffusive regime. 

Under these conditions the random contribution to the 
SO Hamiltonian obtained with Eq.(3) is: 

rr / s , , i ( d« r (p) da r {p) \ 
H T = -a r {p)<7yk x + - [vy—g^ Gx ~Qy~ ) ' 

The spin of a moving electron is, therefore, a subject 
of a randomly time-dependent SO field, with each elec- 
tron probing the realization of the field at the given path 
or different electrons probing the different configurations 
of essentially the same random field. The direction and 
the magnitude of the Zeeman field randomly change in 
time, and therefore, a dephasing of the spin states oc- 
curs. The correlation function of the random field is 
characterized by the Fourier component (a r (0)a r (r)) w 
extended in the range of frequencies lu up to fid ~ l/r^ 
where = zq /v is the time the electron takes to travel a 
distance zq through which the Rashba coupling changes 
appreciably and v is the electron velocity. The fre- 
quency fid is much larger than the spin-flip frequency 
flso ~ ctRk/h with Qso/^d being of the order of 0.1. 
This condition implies that the random contribution to 
the Zeeman field is out of resonance with the spin-flip 
transitions, and, for this reason, the spin-flip transition 
rate is low. Therefore, under these conditions it is neces- 
sary to pass through many domains of the SO coupling 
for S to be destroyed. Analysis similar to that done in 
Refs.mH shows that in the case Qd ^> Hso the ran- 
dom SO coupling mechanism leads to the spin relaxation 
rate of the order of 27r(a 2 )zomax(fc 2 , ZQ 2 )/h 2 v. In the 
system considered here, with k 2 z^ 2 it scales with 
the system parameters as (ur(V = 0)) /nz$k. This trend 
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is seen in Fig. 3 which presents the results of the Monte 
Carlo simulation of the spin dynamics. As one can see 
in Fig. 3, {(a x ) + {a y ))/y/2 is gradually destroyed at the 
electron path L s of the order of 50 /xm, that is about 
5 x 10 3 domains of the SO random coupling. We asso- 
ciate this decay time with the spin-coherence lifetime in 
the system. The role of the randomness is demonstrated 
by the fact that the spin-coherence lifetime depends very 
strongly on the distance between the dopant layer and 
the conducting electrons increasing with the distance Zq. 
Even when the mean SO coupling is the same, the in- 
crease in the distance between the dopant layer and the 
symmetry plane of the quantum well decreases the ran- 
domness and, in turn, the spin relaxation rate, as can 
be seen from a comparison of the S(t) dependences in 
Fig. 3. The noise in the £(£)-dependence corresponds to 
the simulation with a finite number of electrons and must 
be experimentally observable in the SFETs operating in 
this mode. For these reasons, the size of the transistor 
base cannot be very long, and the spins cannot be con- 
trollably delivered at distances much larger than 100 /im 
at best. 

An order of magnitude estimate of the ratio of the re- 
laxation rates due to these two mechanisms, namely, the 
SO coupling randomness T r and the fc 3 — terms contribu- 
tions To gives: 




The small ratio r^/rfc, which strongly depends on the 
mobility, is typically of the order of 10~ 2 , and, there- 
fore, favors the role of the fc 3 Dresselhaus contribution. 
However, at wk < 1, it is compensated by the 7r 4 prefac- 
tor arising due to the size quantization and by the effect 
that due to the fast in-plane angular dependence of the 
fc 3 — originated terms, the efficiency of this mechanism is 
decreased, leading to a large prefactor 16ir in Eq|SJ The 
exact value of this prefactor is model-dependent, being, 
however, of the order of ten at any reasonable model. 
Therefore, even at relatively small fluctuations of the 
Rashba SO coupling (a 2 ) /{a) 2 R ~ 0.1, they will dom- 
inate as the spin relaxation mechanism. We mention 
here that despite the spin-orbit coupling in our model 
depending linearly on the momentum, in the quasibal- 
listic regime the spin relaxation path L s = Hk/mT r , 
where m is electron effective mass, shows only a weak 
electron-momentum dependence. The reason is that 
the random contribution to the spin precession angle 
for passing through one domain, being of the order 
of a T k/h x zo/ihk/m), is momentum-independent, and, 
therefore, the number of the domains contributing into 
L s , depends weakly on k. At the same time, for the Dres- 
selhaus fc 3 — contribution one would expect L s ~ fc~ 5 . 



IV. THE (Oil) QUANTUM WELL 

Another interesting system where the direction of 
the SO field is momentum-independent is the (Oil) 
zincblende quantum well. Here SO coupling is described 
by the Dresselhaus Hamiltonian 

H D = a D k y a z [1 - (fc 2 - 2fc 2 ) /(fc 2 )] , (9) 

with the direction of the SO field being always parallel 
to the z-axis which is perpendicular to the QW plane. 
Here a D = a c (k 2 z )/2. The (fc 2 - 2fc 2 ) /(fc 2 ) term in the 
square brackets is the bulk-originated Dresselhaus term, 
which, contrary to the (001) quantum well, due to the 
symmetry, does not lead to the change in the SO field 
direction with the changes in the momentum. When the 
system is the subject of the regular Rashba Hamiltonian 
Hji = (an) (a x k y — a y k x ) also, the electron spin precesses 
around the axis determined by the direction of the Hd + 

H R field with the rate fi = 2^Ja 2 D k 2 + a 2 R k 2 /h, where 

fc 2 = fc 2 + fc 2 . For the component of the spin initially 
polarized along the z— axis, one obtains 

<J z {t) = ct z + A CT cos fit, (10) 
n, 2 k 2 

° Z ~ n 2 fc 2 -I- it 2 fc 2 ' 

a D K y + a fl «|| 

rv 2 fc 2 

A = a R^\\ 

a 2 D k 2 y + a 2 R k 2 ' 

where A a is the precession amplitude. Randomness in 
Hr causes fluctuations in the magnitude and direction 
of the Zeeman SO field acting on the electron spin, both 
now changing randomly in time. Eventually, this leads 
to the spin relaxation, which can be also understood as 
the a z = 1 to a z — — 1 spin- flip transitions caused by the 
random field. 

We consider the case of the asymmetrically doped 
quantum well for which the mean value of the Rashba 
parameter can be influenced by applied bias. A typical 
evolution of (a z (t)) averaged over an ensemble of elec- 
trons for the parameters of spin-orbit coupling typical 
for Ino.5Gao.5As quantum wells is presented in Fig. 4. 
We consider two cases: a partially and the fully compen- 
sated Rashba terms. In the former case the relaxation 
occurs as both the amplitude of the oscillations and the 
mean value averaged over the ensemble of electrons tend 
to zero. As one can see in the Fig. 4, the decay time of 
(&z{t)) m the Ino.5Gao.5As structures is of the order of 10 
ps, similar to the case of the ao = (or(V)) transistor. 

In "pure" GaAs structures, where a$o is an order of 
magnitude smaller than in the Ino.5Gao.5As alloys, one 
would expect the spin relaxation time due to the ran- 
domness of the doping which roughly scales as otgQ, two 
orders of magnitude longer, that is of the order of 10 3 ps. 
A long electron spin coherence time T s of this order of 
magnitude has indeed been observed experimentally in 
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Refs.|27j for the (Oil) symmetrically doped GaAs QWs 
of relatively low mobility with fj, ~ 10 3 cm 2 /(Vs). It 
is instructive to estimate the corresponding distance L a 
at which the spins of diffusively propagating electrons 
can be delivered for this time in these low-mobility sam- 
ples. With this mobility, the momentum relaxation time 
t ~ 0.1 ps leads to L s ~ v^t) 1 ' 2 , where v ~ 10 7 cm/s 
is the speed of optically injected electrons investigated 
in Ref. 27]. Here (T^r) 1 / 2 - 10 ps, and therefore, L s 
is of the order of a few /im. Therefore, due to the low 
mobility these QWs cannot be used for the simultane- 
ous long-distance charge and spin transfer necessary for 
the spintronics applications despite a long spin coherence 
time there. 



V. CONCLUSIONS 

We have investigated the spin relaxation rate arising 
due to a random SO coupling in a possible realizations 
of a SFET with the balanced Rashba and Dresselhaus 
spin-orbit coupling and the (Oil) zincblende quantum 
well. These systems, where the direction of the effective 



spin-orbit Zeeman field is expected to be independent on 
the electron momentum and, therefore, the Dyakonov- 
Perel mechanism of spin relaxation must be suppressed, 
are considered as promising elements for the spintronics 
applications. At the same time, due to the randomness 
of the Rashba contribution to the spin-orbit coupling, 
causing an additional spin relaxation, in both cases the 
relaxation time cannot exceed considerably 10 ps thus 
limiting the size of the base of this transistor of the order 
of 10 /xm, which weakly depends on the electron mo- 
mentum, and, therefore, restricts the possibility of their 
experimental realization. Recent focus in spintronics is 
related to Si/Ge two-dimensional structures, where the 
spin-orbit coupling is weak 22 - 28 ! 2 ^ 3 ^. Within these sys- 
tems the randomness of the SO coupling, which is re- 
lated either to the randomness in the dopant distribu- 
tion as considered here or to random bonds on the Si/Ge 
interfaces^ will cause spin relaxation and limit their po- 
tential applications as well. 
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FIG. 1: A plot of a typical realization of the random Rashba SO coupling (in arbitrary units) arising due to random variations 
in the dopant concentration, zo = 20 nm, and the size of the template is 200 x 200 nm. The quantum point contacts are 
marked by black triangles, thin lines show the electron paths, and black arrows show the direction of the electron spins. 
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FIG. 2: Typical realization of the random Rashba SO coupling for different distances between the dopant layer and the 
symmetry plane of the QW as marked near the plots, n = 2.5 x 10 11 cm~ 2 . 
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FIG. 3: Dependence of {a x +a y } j\[2 on the electron path L, (an{V = 0)} = 6.0 x 1CF 10 eVcm. The simulation includes N = 64 
electrons. It is assumed that the regular part of the Rashba term is reduced to 0.6 of its initial value ((an(V)} = 0.6(au(V = 0)}) 
by the applied bias. The wavevector of electron k = 0.5 x 10 6 cm -1 , and the concentration of dopant ions is n — 5 x 10 11 cm -2 . 
The inset shows spin relaxation aXn — 2.5 x 10 11 cm~ 2 , zo = 8 nm. 




FIG. 4: Time dependence of (<j z ) for the (Oil) quantum well. The distance between the QW symmetry plane and dopant layer 
is 10 nm. The Dresselhaus SO coupling cid = 2.0 x 10~ 10 eVcm corresponds to an Ino.5Gao.5As with a c = 25 eVA 3 QW with 
w = 80 A. The wavevector of electron k y = 2 x 10 6 cm -1 , k x = 0, n = 2.5 x 10 11 cm -2 , and (a R {V = 0)} = 3 x 10 -10 eVcm. 
Dashed line represents {aii(V)) = (an(V = 0))/2, solid line represents (an(V)} = (full compensation). 



